One-Dimensional Multi-Band Correlated Conductors and Anderson Impurity Physics 
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A single Anderson impurity model recently predicted, through its unstable fixed point, the phase 
diagram of a two band model correlated conductor, well confirmed by Dynamical Mean Field Theory 
in infinite dimensions. We study here the one dimensional version of the same model and extract 
its phase diagram in this opposite limit of reduced dimensionality. As expected for one dimension, 
the Mott metal-insulator transition at half filling is replaced by a dimerized insulator-undimerized 
Mott insulator transition, while away from half filling the strongly correlated superconductivity for 
inverted Hund's rule exchange in infinite dimensions is replaced by dominant pairing fluctuations. 
Many other aspects of the one dimensional system, in particular the field theories and their sym- 
metries are remarkably the same as those of the Anderson impurity, whose importance appears 
enhanced. 

PACS numbers: 71.30+h, 71.10.Pm, 71.10.Fd 



Very substantial progress in our understanding of the 
Mott metal-insulator transition (MIT) have been made 
thanks to the so-called Dynamical Mean Field Theory 
(DMFT)P|, a quantum analogue of the classical mean- 
field theory which treats time correlations and is exact 
in infinite dimensions (oo-D). In DMFT, the approach 
to the MIT from the metal phase is accompanied by a 
net separation of energy scales between well pre-formed 
high energy Hubbard bands - images of the excitations in 
the nearby Mott insulator - and the lingering low-energy 
itinerant quasiparticles. This separation is in fact al- 
ready contained in the isolated Anderson impurity model 
(AIM), where most of the spectral weight is concentrated 
in the high-energy sub-bands and only a small fraction 
- describing quasiparticles promoted into the conduction 
screening-bath - remains close to the chemical poten- 
tial. This is of course unsurprising since in oo-D DMFT 
maps the lattice model of interacting electrons onto an 
AIM supplemented by a self-consistency condition[l|. Ir- 
respective of whether this mapping is merely a trick to 
solve the lattice model in oo-D or whether it hides per- 
haps a more fundamental aspect of the physics close to a 
MIT, this does suggest that some of the strongly corre- 
lated lattice properties could be directly inferred by the 
AIM itself, even without self-consistency. This route was 
recently explored to anticipate the anomalous properties 
near the MIT of a two-band Hubbard conductor on the 
basis of the phase diagram and in particular of the unsta- 
ble fixed point of a two-orbital AIM|2(. All the predicted 
properties, including strongly correlated superconductiv- 
ity near the Mott transition Q were later confirmed by 
full DMFT0]. Despite that success, it would still seem 
hazardous to suggest that the properties of a single AIM 
have generally anything to do with the actual behavior of 
the model lattice conductor away from oo-D, least of all 
in the opposite extreme of one dimension (ID). We show 
in this Letter that, apart from some obvious differences 
related to dimensionality, the phase diagram of the model 
does not change significantly in ID, where therefore the 



AIM physics appears to remain significant. 

We consider the two-band Hamiltonian near half-filling 
and in ID 

H = t £ £ (4, aCT c i+liOCT + H.c) 2) 2 
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. , aa ) creates(annihilates) an electron at site 
/, in orbital a — 1,2 with spin a =T,|, and hi — 



J2aa c l,aa c i,aa * s * ne electron density at site i, and U ^> 
|J| is an on-site Coulomb repulsion. This Hamiltonian 
has recently been discussed as relevant to the novel doped 
metal-phthalocyanine conductors. For J = Eq. Q 
describes an SU(4) Hubbard model, analysed e.g. in 
Ref . [|| . A finite value of J lowers the symmetry down to 
U(l)xSU(2)x (U(l) x Z 2 ), where U(l) refers to charge, 
SU(2) to spin and (U(l) x Z 2 ) to the flavour (orbital) 
sector. We stress here that the single AIM shares identi- 
cally this same symmetry, a point which we will return to 
further down. Two electrons on the same site can form 
either a spin triplet, with energy J/4, an inter-orbital 
singlet, with energy —3 J/4, or two intra-orbital singlets 
with energy — J/4. Therefore J < favors the spin triplet 
while J > the inter-orbital singlet. Actually J > 0, 
("inverted Hund's rule exchange"), provides a pairing 
mechanism in the Cooper channel cj^c^, + c\ 2 jcj im- 
pairing is impeded by the repulsion U, so that the bare 
scattering amplitude in the inter-chain singlet channel, 
A = U — J/2, is attractive only in the unrealistic case 
of J > 2U > 0, apparently excluding superconductiv- 
ity despite the pairing mechanism provided by J > 0. 
As shown in DMFT Q, this naive expection is actually 
wrong, at least in oo-D. A superconducting pocket ap- 
pears near the half-filled Mott insulator, U ~ t ^> J, 
moreover with a hugely enhanced superconducting gap 
with respect to the U = BCS gap value. This surpris- 
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ing result had in fact been foreshadowed by the single 
AIM study whose phase diagram displays an unstable 
fixed point at < J* ~ Tr-, where Tk is the Kondo tem- 
perature, which separates a Kondo screened phase for 
J < J* from an unscreened phase for J > J» . Translated 
into DMFT, the AIM Kondo temperature becomes the 
quasiparticle coherent bandwidth, vanishing at the MIT. 
This implies that the AIM onto which the metallic lattice 
model maps in oo-D must necessarily cross the unstable 
fixed point J ~ Tk before the MIT. The speculation Q 
that the lattice model would respond to the local insta- 
bility by spontaneously developing a bulk order parame- 
ter in the inter-orbital singlet Cooper channel was fully 
confirmed by DMFT [3]. 

We turn now to study the same model in ID. As 
usual it is convenient here to represent Q within 
bosonization. The Fermi fields around the right (R), 
+kp, and left (L), —hp, Fermi points are expressed as 

i>R(L),aa{x) ~ exp [±lk F X - 1 y^TT (6 aCT ( x) =F ao (x) )] , 

where ^aa{x) and d x Q a a{x) are conjugate 
Bose fields. We introduce the linear Com- 

mit +$ U +$2 T +$2l)/2, 



($i T - $ u +$2T - *2i)/2, $ 



binations 
$ s 

($ 1T + $ u - $ 2T - $ 2 j.)/2, and $ s/ 
(<&ij — — + $24.) /2, which describe respec- 
tively the total charge, the total spin, the relative 
charge and the relative spin density fluctuations. In this 
representation the interaction involves only bilinears of 
cos \/4jr<b n and cos ViirQn, n — c,s, /, sf, which in turn 
can be expressed as 

COsV47T$„ = -i— {£R, n £L,n + (R,n(L,n) , (2) 
COsV47T0„ = i— (%R, n €L,n - (R,n(L,n) > (3) 

where £,R(L),n and Cr(l),u are Majorana fermions and a 
is a cutoff distance. These fermions can be used to intro- 
duce eight two-dimensional classical Ising models, each 
one in principle characterized by a mass to ~ (T — T c ), 
to < and m > meaning ordered and disordered 
phases, and to = the critical point. In this way the U(l) 
charge sector is represented by two identical Ising models 
with mass m c , the doublet (£r iC , £ i)C ) and (Cr, c , Cl,c); the 
SU(2) spin sector by three identical Ising models (mass 
m s ), the triplet {£r, s ,£l,s), {Cr,s,Cl,s) and £z, a /); 
the U(l) flavour sector by two identical Ising copies (mass 
to/), the doublet (£r,/,£l,/) and (CrjXlj); and fi- 
nally the remaining flavour Z2 by a single Ising model 
(mass mo), the singlet (Cr,s/, &,«/)• Hi Without inter- 
action all Ising models are critical. The interaction in- 
duces marginally relevant couplings between them which 
might spontaneously generate finite masses. Indeed, by 
a fermion two-loop renormalization group (RG) analysis, 
we find that the Hamiltonian generally flows to strong 
coupling fixed points which allow a simple mean-field de- 
scription in terms of finite average values of ((,R t n^L,n) 
and (Cii,nCi,n)) n = c i s if> s fi which preserve all con- 
tinuous symmetries. Within this same description the 



phase diagram of (JIJ turns out to be characterized by 
simply identifying the relative signs of the masses rrzj, 
i = c,s, /, 0, while the overall sign has a physical meanin 
only in the spontaneoulsy dimcrized phase, see below 

We start from half filling, (hi) = 2, and analyse the 
phase diagram for increasing U ft keeping for simplicity 
a fixed ratio U/\J\ 3> 1. At weak coupling, U/t <§; 1, the 
Hamiltonian flows under RG to an SO (8) Gross-Neveu 
model, which describes a spontaneously dimerized in- 
sulator with gaps in the whole excitation spectrum. It 
is known that the SU(4) Hubbard model, i.e. JTJ with 
J = 0, dimerizes at half- filling [3. ITo|. and that a small 
I J\ <C U cannot destabilize this gapped phase. This 
phase is characterized by all masses having the same sign, 
with the overall sign reflecting the broken translational 
symmetry. Eventually, though, this dimerized phase can- 
not survive indefinitely for large U/t. When U ^> t, two 
electrons localize at each site in a configuration optimiz- 
ing the on-site exchange J. In particular for J < 0, con- 
ventional Hund's rules, the model effectively reduces to 
a spin S=l Hcisenberg chain, still gapped everywhere in 
the spectrum but not dimerized|l2|. For inverted Hund's 
rules, J > , the singlet configuration 



C l, lt C l,2J. + C i,2T C i,lJ. 



|0>, 
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is favored and the ground state is akin to a collection of 
local singlets, a kind of local valence-bond Mott insula- 
tor, still gapped but not dimerized. We conclude that 
upon increasing U/t the dimerization must disappear for 
either sign of J. In fact, whereas at weak coupling at 
for J = the dimerization-induced gaps in the spin and 
flavour sectors follow the BCS-like behavior of the charge 
gap, the latter continues to increase monotonically as U 
increases while the former gaps reach a maximum, ap- 
proximately when U ~ 5£, and then start dropping as 
t 2 /U for U » This decoupling of charge from 

spin and orbital modes is a ID remnant of the MIT and 
seems quite sharp [(J. The weak exchange \J\/U <C 1 , 
irrelevant in the weak coupling dimerized phase, eventu- 
ally turns in strong coupling to a relevant perturbation 
able to suppress dimerization. Thus one might expect 
that this could occur only for very large U, when the 
spin gap induced by dimerization becomes small of order 
I J\. However it cannot be excluded that the demise of 
dimerization could even take place for smaller U , say for 
U < 5t. The two-loop RG equations moreover suggest 
for J < a c = 3/2 spin-SU(2)2 critical point where 
the triplet mass m s changes sign 13], signaling a transi- 
tion from the dimerized insulator to the Haldane spin-1 
chain Mott insulator. For J > the transition is in- 
stead predicted to occur through a c = 1/2 Ising critical 
point where the singlet mass too crosses zero, signaling 
the transition from a dimerized to the valence-bond Mott 
insulator, see Fig. ^ 

Let us introduce doping, moving away from half filling. 
Dropping the Umklapp terms from the weak coupling 
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Valence-bond MI Dimerized MI Haldane MI 



FIG. 1: Phase diagram of model Q near half- filling as a func- 
tion of doping \n - 2| and of U/t at fixed | J\/U < 1 for J < 
(right panel) and J > (left panel) . At half-filling the system 
is always insulating and displays by increasing U/t an Ising 
transition (c = 1/2 flavour) from a Dimerized Mott insula- 
tor (MI) to a Valence-bond MI for J > 0; and a c = 3/2 
transition to a S=l Haldane MI for J < where the spin is 
gapless (c — 3/2 spin). All the insulating phases evolve upon 
doping into metal with gaps in all noncharge sectors. For the 
doped Haldane MI the additional label (flavour) or (spin) in- 
dicates that the phase can be viewed as the natural evolution 
away from half-filling of a Haldane chain built of either spin 
or flavour triplets. The transition lines bewteen each different 
metallic phase are identified by the central charge c and by 
the sector involved, spin or flavour. 



RG equations, the interaction flows, for either sign of J, 
towards a fixed point 

H in t -> ~g* J dx A f {x) AO), (5) 

where 

(6) 

which is a spin-singlet, flavour-singlet but space-odd pair- 
ing operator ^3- This fixed point interaction has a 
dynamically enlarged SU(4) symmetry [{j, ITU . unlike the 
original model and realizes a doped Haldane chain, 
the two 1/2-spin constituents forming sing let bonds, one 
to its right and the other to its left.[15| Moving from 
J < to J > interchanges the spin sector and the 
flavour sector, which includes the doublet and the sin- 
glet forming a degenerate triplet, again a dynamically en- 
larged SU(2)2 flavour symmetry. The Ising masses satisfy 
m c = 0, nif Too > but m s m/ < 0. The pairing correla- 
tion function (A(x)A+(0)) ~ {l/x) 1/2K " where K c is the 
Luttingcr liquid exponent of the gapless U(l) charge sec- 
tor. There are also power-law decaying correlation 
functions with exponent 2K C which involve the density- 
wave operators exp (±i4 kp x ± iv'^'&c) cos \/47r$ s (/) 
and exp (±iA k F x ± iViTT^ c ) cos\/47r9 s /. If K c > 1/2 
the superconducting fluctuations dominate over the Akp 
density-wave ones, and the opposite for K c < 1/2. Since 
the model has an insulating phase at quarter filling [l6j. 
by standard arguments[l7j we expect K c > 1/4. Hence 



the pairing susceptibility in channel © always diverges 
faster than for free fermions. 

Revealing as it is, this weak coupling analysis is not 
fully satisfying as it implies an unphysically abrupt 
change of sign of m s at the slightest density deviation 
from half-filling. A better approach near half filling may 
be a two-cutoff RG scheme, namely running at first the 
RG as if for half filling until reaching an energy scale of 
the order of the chemical potential shift, and only at this 
point dropping the Umklapp terms. Doing this we find 
that the doped dimerized Mott insulator turns to a metal- 
lic phase, m c = 0, all other masses retaining their sign. 
Here the dimer order parameter is zero but its correlation 
function decays slowly with a power-law exponent K c /2. 
This agrees with a similar analysis by Boulat^O] and 
leads us to propose the phase diagram of Fig. ^where the 
doped dimerized insulator transform into the doped Hal- 
dane phase by a c = 3/2 critical line^j| which involves 
the spin or the flavour sector for J < and J > 0, respec- 
tively. In addition we expect, by similar arguments, that 
the valence bond Mott insulator at large U/t with J > 
transforms upon doping to a metal phase, with m c = 0, 
m s mf > but momj < 0. This phase is identified by 
a fixed point interaction of exactly the same form as (JSJ) 
with the pairing operator corresponding to the singlet 
configuration namely 

(7) 

which is still a spin-singlet but it is now no longer in- 
variant under the flavour SU(2). The pairing correlation 
function in this channel still decays with exponent l/2K Cl 
and again there are competing Akp density wave fluctu- 
ations with exponent 2K C . Therefore we argue that the 
Ising critical point at half-filling for J > and large U/t 
extends to a critical c = 1/2 line away from half-filling, 
as in Fig. 2] This line should merge into the c = 3/2 
flavor SU(2)2 critical line, out of which a c = 1 U(l) crit- 
ical line must emerge. Along this line the flavor doublet 
mass, to/, changes sign. This scenario is compatible with 
our expectation far away from half-filling. Indeed, if we 
keep the ratio J/U > fixed and increase U/t we arrive 
at a situation where J>f. Here, whenever two electrons 
occupy the same site, they are forced into the singlet con- 
figuration Q. This constraint can be implemented by a 
projector lea ding to a model which was numerically anal- 
ysed in Ref. Jlj| close to quarter filling. The numerical 
results are compatible with the existence of the doped 
valence-bond phase. Since the weak coupling phase is in- 
stead the doped Haldane chain, we conclude that there 
is a transition by increasing U/t at fixed J/U > whose 
criticality belongs to the U(l) universality class. Out 
of all these arguments we finally draw the quantitative 
phase diagram of Fig. ^ 

We note the presence throughout the phase diagram 
(except of course at half filling) of a singular spin-singlet 
Cooper pairing susceptibility either in channel JSJ, within 
the doped Haldane insulator, or in channel Q, within the 
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doped valence-bond insulator. The latter, as discussed at 
the beginning, is most unexpected, emerging only for suf- 
ficiently strong repulsion U. We note that, if < J <C t 
were fixed, we still would expect at U = singular pairing 
susceptibilities below an exponentially small energy scale 
of the order of the spin and flavour gaps. Remarkably, we 
find that the role of a large U is to raise this energy scale 
by increasing the magnitude of the gaps, which in turns 
implies that U effectively enhances pairing fluctuations. 

Finally we consider the competition between Cooper 
pairing and Akp density- wave fluctuations in the doped 
valence bond and in the Haldane phases as we approach 
the Mott insulator at half filling. The insulating behavior 
is driven at weak coupling by 4/cF-Umklapp terms of the 
form cos\/47r^c cos \/~4jt<& s (cos\/47r$/) (cos v4tt0s/) ■ 
Since the spin and flavour sectors remain gapped away 
from half filling, one is tempted to conclude that K c — > 1 
on approaching half filling[l7|- However, both Haldane 
and valence bond Mott insulators require a large U/t, and 
the above weak-coupling argument might not be correct. 
One may argue that the decoupling of charge from the 
other sectors as U/t increases may be modelled by adding 
an 8k p -Umklapp cos -\/167r$ c @- If alone this term would 
rather suggest a MIT at half filling with K c — > 1/4. How- 
ever it was found numerically that in the half-filled SU(4) 
Hubbard model the charge-2 Majorana fermions, £r(l),c 
and Cr(l),cj remain coherent excitations even at large 
U, even though their energy is only slightly below the 
two-particle continuum, merging into the latter only for 
U/t — ► oo. This suggests that K c tends always to 1 as 
the density approaches half filling, even though a larger 
U implies a sharper crossover from K c ~ 1/4 to K c = 1 . 
Hence we conclude that sufficiently close to the MIT the 
Cooper channel pairing dominates over Akp fluctuations. 

We can finally return to our original motivation and 
discuss differences/analogies with the phase diagram of 
the same model suggested by the AIM and verified by 
DMFT in infinite dimensions for J > 0. At half fill- 



ing we always find in ID an insulator because of perfect 
nesting, generally absent in higher dimensions. Hence 
the MIT of oo-D is replaced in ID by an Ising transition 
between two insulators, the band-like dimerized insulator 
(driven by nesting) and the strong-coupling undimerized 
valence-bond Mott insulator. Upon doping, we still ex- 
pect by increasing U /t an Ising transition from a doped 
dimerized insulator into a doped valence-bond insulator, 
the latter characterized by a singular pairing susceptibil- 
ity. This is the exact analog of what was found in infinite 
dimensions. Even more surprising is the role of the sin- 
glet Ising sector ((fl,s/, Ci,s/) which becomes critical at 
the transition. As previously mentioned, the behavior of 
the lattice model in oo-D is controlled near the MIT by 
the unstable fixed point of the AIM onto which the lat- 
tice model maps. In turn this unstable fixed point can 
be interpreted as the free boundary condition fixed point 
separating the two fixed boundary conditions just in the 
same Ising sector, and that is the natural generalization 
to a boundary problem of the Ising critical point which 
we uncovered in ID. It is now clear that both the ID and 
the single impurity analysis make use of the same field 
theoretical scheme, the non-abelian bosonization, where 
the Ising sector emerges naturally from embedding the 
flavour SU(2) 2 into U(l). Suggestively, this appears to 
be the reason why the behavior in one dimension and in 
infinite dimensions are so similar. At this stage we cannot 
say whether this similarity is a mere accident or not. It 
certainly does encourage the speculation that the physics 
of a single AIM may play a more fundamental role in the 
description of strongly correlated metals in any dimen- 
sions, at least at intermediate energy/temperature scales 
before full bulk coherence settles in. 
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